Correlated and Decorrelated Positional and 
Orientational Order in the Nucleosomal Core 
Particle Mesophases 

V. Lorman \ R. Podgornik ^'^'^ and B. Zeks ^'^ 
^ Laboratoire de Physique Mathematique et Theorique, 
Universite Montpellier II, F-34095 Montpellier, France 
^ Department of Physics, Faculty of Mathematics and Physics, 
University of Ljubljana, SI- 1000 Ljubljana, Slovenia 
^ Department of Theoretical Physics, J.Stefan Institute, 
SI-1000 Ljubljana, Slovenia 
4 LPSB/NICHHD, Bid. 9 Rm. lE-116, National Institutes of Health, 

Bethesda, MD 20892-5626 
^ Institute of Biophysics, Medical Faculty, 
University of Ljubljana, SI-1000 Ljubljana, Slovenia 

February 2, 2008 



Abstract 

We investigate the orientational order of transverse polarization 
vectors of long columns of nucleosomal core particles and their cou- 
pling to positional order in high density mesophases discovered re- 
cently. Inhomogeneous polar ordering of these columns precipitates 
crystalization of the 2D sections with different orientations of the 
transverse polarization vector on each column in the unit cell. We 
propose possible scenarios for going from the 2D hexagonal phase into 
distorted lamellar and related phases observed experimentally. 



The problem of efficient genome compaction in the evolution of eucary- 
otic organisms was solved on the smallest scale via the formation of DNA 



- histone protein complexes, that make a fundamental unit of chromatin: 
the nucleosome PP. A nucleosome is composed of a 147 bp long DNA frag- 
ment, wrapped 1.75 times around the histone octamer core composed of four 
different histone proteins and a DNA linker that connects two consecutive 
nucleosomes Digestion of the linker DNA between different nucleosomes 
along the eucaryotic genome gives rise to the nucleosome core particle (NCP) . 
NCP has an approximately cylindrical shape of height and radius ~ 55 A, 
and a structural charge of ~ -165. This complex is stable in aqueous solutions 
from mM to 750 mM monovalent salt ionic strength 

At high enough concentration the solution of NCPs ceases to be isotropic. 
Competing, as yet poorly understood [H |H] , microscopic interactions between 
NCPs in concentrated solutions give rise to a slew of liquid crystalline phases 
of fascinating complexity [S' . At low salt concentrations and sufficiently high 
osmotic pressure NCPs aggregate into columns of almost ideal cylindrical 
shape which contain dozens of individual particles. These columns then 
aggregate to form a lamellar phase [7j. At higher salt concentrations this 
phase gives way to a 2D hexagonal columnar liquid crystalline phase leading 
at still higher salt concentrations to a 3D hexagonal crystal [,8J. 

The high asymmetry of the NCP particle, with a pronounced dyadic axis 
marking the entry-exit point of the wrapped nucleosomal DNA, allows one to 
actually observe the local orientational order of the NCPs, described by its 
in-plane dyadic axis vector, conventionally called transverse polarization [0]. 
In the lamellar phase the dyadic axes of individual NCPs are strongly corre- 
lated along the column leading to a non-zero average transverse polarization 
vector P of a column. Cryomicroscopy studies [7] give strong evidence that 
the dyadic axes, and consequently the vectors P, of columns in the two layers 
of a lamella are in an antiparallel direction. In addition, the columns in the 
two layers are shifted to a small but constant distance in such a way that the 
antiparallel dyadic axes form a periodic polar vector field (Fig. 1). SAXS 
refine this picture [H] by showing that there are strong correlations between 
the NCP columns within a single lamella but different lamellae appear to be 
completely decorrelated (Fig. 1). In this sense the lamellar phase is decorre- 
lated at a macroscopic scale and corresponds to simple smectic ordering. 2D 
hexagonal phase is characterized, in turn, by a positional order in the plane 
perpendicular to the columns, but dyadic axes of NCPs are orient at ionally 
disordered along the columns (Fig. 1). The corresponding average polar 
vector of the column is zero. 

In what follows we will investigate the nature of positional and orienta- 
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tional order in phases that are induced by a condensation of the periodic 
antiparallel transverse polarization field modes. Basing our assumptions on 
experiments we postulate that the transverse polarization order in a lamella 
described with P(p), where p = {x,y), is periodic and highly anisotropic. 
We will show that the transition between the 2D colloidal crystal with dis- 
ordered dyadic axes and the smectic-like lamellar phase with the in-plane 
periodic polar order can be understood by introducing a common subphase 
substructure for both phases. The common subphase is a 2D colloidal crystal 
with lamellar structure and with the in-plane polar order of individual lamel- 
lae similar to that in the actual lamellar phase. However, the correlations 
between the positions of columns and the orientations of their dyadic axes 
from one lamella to another make this phase crystalline and different with 
respect to the smectic-like structure. The transformation from the hexagonal 
colloidal crystal to the smectic lamellar phase can take place as a sequence of 
two second-order transitions through the intermediate crystal lamellar phase, 
or as a direct first-order transition. Order parameters of both transforma- 
tions can be easily expressed in terms of the periodic transverse polarization 
fields. 

If the transverse polarization vector field condensation starts from a 2D 
hexagonal phase then in general one would have several different choices for 
the direction of the wave vector k with respect to the symmetry planes of the 
2D hexagonal phase. In the hexagonal coordinate system the wave vector 
of the condensation into the crystal lamellar phase has to be described in 
Kovalev's notation ^U] by ki2 = |bi, where bi is one of the reciprocal 
basis vectors. The irreducible star of the ki2 vector has three arms given 
by ki = k^^j = |bi, k2 = kS? = ib2 and kg = kj? = i(-bi + b2). 
The basis functions of a polar vector field with this periodicity are given by 
^1 = {2e, + ey) e^i'i ^ ^2 = (e. + 2e,) e'''^ '^ and ^3 = (-e^ + e^) e^'^^'^ 
The physical realisation of the order parameter is given by the symmetry 
coordinates composed of the column's polar vectors Pj. Denoting by i the 
index of the column, in such a way that its (x, y) coordinates are i = 1 — > 
(0,0), i = 2 — > (1,0), i = 3 — > (0,1) and i = 4 — > (1,1) in the 2D 
hexagonal unit cell and the 2D polar vector components by {Px,Py), one 
gets the following components of the order parameter 

Vl = Pxl~Px2+Px3~Px4r\-\ {.Pyl—Py2+Py3~Pyi) 
^2 = \ {Pxl+Px2—Px3—Px4)+Pyl+Py2—Py3—Py4: 

^3 = 1 ^xl+Px2+Px3~Px4+Pyl—Py2~Py3+Py4:) ■ (l) 
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The Landau free energy of this system is now given via the orientation prob- 
abihty density function p(x, y, z) = pq + Yl^i=i Vii'ii^^ the form 

•^ = y" f {Vi'^iix,y,z)) dxdydz = F{r]i). (2) 

For the sake of clarity, we have omitted in Eq. |21the terms stemming from the 
large-scale elastic deformations of this order. Consequently, the free energy 
of the system can be considered as a simple function of the order parameter 
components. 

Let us also stress an evident duality between the polar vector carried 
by a NCP column and the polar vector of its displacement in the solvent 
which has important consequencies for the structure of its ordered phases. 
Since the displacements Uj of NCP columns and the polar vectors Pj of the 
dyadic axes span the same irreducible representation of the symmetry group 
of the 2D hexagonal phase the components of the order parameter can be 
expressed in the same form in terms of the columns' displacements {u^i, Uyi). 
Minimization of the free energy dependent on both Pj and Uj shows that 
they are proportional. This fact permits one to obtain also the positional 
order and thus also the final structures of the phases. 

To establish the exaustive list of the phases induced by the considered 
condensation mechanism we use the fact that F{rii) has to be an invariant 
function of the irreducible representation spanned by the periodic transverse 
polarization vector field. Since the basis invariants of this represenatation 
are Ii = rjf + + r]'^, I2 = '7i'72 +^2^3 + ^3^1 I3 = r(lrjlr(l, the free energy 
Eq. 121 has to have the form JF = F (Ji, J2, /s). The minimization of this free 
energy w.r.t rji gives the equation of state. There are six different solutions 
of this equation of state which we list below in terms of r/j, the unit cell 
multiplication in the ordered state with respect to the 2D hexagonal phase 
i.e. Vord/Vhex and their space groups. Let us show that the first solution, 
i.e. rji 7^ 0, ri2;i = describes the crystal lamellar state. From Eq. ^ it 
follows that in this case Pi = P3 and P2 = P4, or in other words, that the 
orientations of dyadic axes carried by columns 1 and 3 or columns 2 and 4 
are parellel and have the same degree of correlation along the column. Thus, 
due to the duality, the displacements of the columns 1 and 3 as well as 2 and 
4 are also identical. This means that the periodic transverse polarization 
field induces a lamellar structure, where a single lamella is defined by the 
columns 1,3,1,3... for the top part and 2,4,2,4... for the bottom part of a 
single lamella. The order parameter in this crystalline lamellar phase induced 
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Table 1: Different solutions of the equation of state. 

by the periodic transverse polarization vector ordering is given by 

771 = (2P,i + P,i) - (2P,2 + Pv2). (3) 

The simplest phenomenological Landau free energy expansion of an isolated 
transition from a hexagonal 2D crystal rji = (0, 0, 0) to a lamellar 2D crystal 
with Tji = {fj, 0, 0), with basis invariants Ji = 77^, I2 = rj'^ and I3 = rf can be 
written as 

= aJi + 61/2 + ... . (4) 

This free energy describes a second order transition between a hexagonal and 
a lamellar crystalline phases, driven trivially by the parameter oi, that can 
depend on various solution conditions, e.g. the salt concentration. The same 
periodic transverse polarization vector field can induce also more complicated 
transitions. Assuming now a general form of the order parameter we can start 
with the following free energy ansatz 

J'liVi) ^ aJi + a2l^ + hh + (5) 

The phase diagram in this case is presented in Fig. 2. One can observe that 
in addition to the lamellar phase, stable for bi > 0, the diagram shows the 
(?7,?7,0) phase, to be discussed below, stable for bi < 0. The growth of the 
complexity of the phase diagram as more non-linear terms are included into 
the Landau free energy is clearly apparent. E.g. if one starts from the free 
energy 

^2{Vi) = ^liVi) + asif + cih + ^12/1/2 + ...., (6) 

this leads to the form of the phase diagram presented on Fig. 2, where an 
additional phase, this time (771,772,0), intercalates between the (77,77,0) and 
the (77, 0, 0) phases. 



5 



Figure 1: Crystalline lamellar phase (c) as a common subphase phase of the 
2D hexagonal columnar (a) and the decorrelated smectic lamellar phase (b). 
Effect of chirality on the S-axes of the NCPs (d). The red circles and black 
squares denote the front (F) and the back (B) direction of the NCP dyadic 
axis 6j (lower right). 

The structure of other phases listed in Table [TJ induced by the periodic 
ordering of the transverse polarization vector field obtained in highly non- 
linear free energy models can be quite complicated. We have already analysed 
the lamellar phase, beside this we also have: 

1. {r],r],ri). The corresponding vector field is given in this phase by the 
vectors of the different columns in the form : 

Pi = (P,,2P,) P2 = (0,0); 

P3 = (P., -Py) P4 = (-2P., -Py) (7) 

The vectors of orientation (and, consequently, the displacements) of the 
columns are expressed in terms of one polar vector [P^^Py). Note that, be- 
cause of the specific interaction between the columns in this phase, column 
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2 is frustrated : its average polar vector is zero. Consequently, its displace- 
ment with respect to the position in the 2D hexagonal phase is zero and the 
dyadic axes of individual NCPs are disordered along this column. All three 
components of the order parameter are equal to the following combination 
of Pj vectors : rj = P^i + 2Pyi + P^s — Pys — IP^^ — Py4. Analyzing the 
displacement field in this phase, see Fig. 2, it is obvious that the columns of 
a hexagonal phase are now displaced to make an inverted hexagonal phase, 
observed also in experiments whose vertices are given by the columns 3', 
2, 1', 2, 4' and 2. 

2. {ri,ri,0). The corresponding vector displacement fields in this phase are 
given by 



Both non-zero components of the order parameter in this phase are equal 



to r/ = |(P,i + Pyi) + i(-P,2 + Py2) + \{P.z - Pys) ' U^-^ + ^v^) ' ^he 



columns of a hexagonal phase undergoing the orientational and displacement 
fields, see Fig. 2, form now an orthorhombic phase with a large unit cell. It 
can be understood as a succession of two series of columns undergoing ID 
displacement waves : in the first series the wave is transeverse and in the 
second one it is longitudinal. Note that the point symmetry of this structure 
shows only slight deviation from the hexagonal symmetry. 
3. Along the same line the vector fields and the structures of the low- 
symmetry phases (?7i, ''725 0), (r^i, r^i, ?72) and ■^2, can be obtained. 

We have shown how a 2D ordering of a transverse vector polarization 
field can induce a transition into a lamellar crystalline phase. The experi- 
mentally observed lamellar phase is however not of this type and is basically 
equivalent to a smectic phase. The difference is subtle and referrs to the 
fact that in SAXS in addition to the smectic peaks only reflections with two 
non-zero indices from the NCPs correlated along the column and within a 
single lamella are present. No reflection with three non-zero indices which 
can evidence the correlation between different lamellae is observed. The 
lamellae are thus decorrelated and the 3D long range order of the trans- 
verse polarization vector is destroyed. This means that the macroscopic 
symmetry of this phase is higher than that of the crystal lamellar phase with 
1i = {fli 0) 0) and equivalent to the symmetry of a usual lamellar smectic 
phase {G = D^oh ® T^). One can speculate that this decorrelation is related 



P3=^ 
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to the difference in correlation lengths of microscopic interactions which in- 
sure the stability of NCP mesophases on the one hand and their colloidal 
character on the other hand. Indeed, the displacements of columns in a 2D 
colloidal crystal are much greater than typical values of atomic shifts in a 
close-packed hexagonal crystal. Resulting distance di between the lamellae 
in our case can become greater than the correlation length of the interaction 
between the NCPs but still smaller then the correlation length of the smectic 
density wave ^p{z)- Since the symmetry of the crystal lamellar phase is lower 




(i^.o.oj (n,ii,0) tn.n-n) 




Figure 2: Displacement fields for the transition from the hexagonal to lamel- 
lar 0, 0) (columns forming a single lamella are enclosed within a gray rect- 
angle), orthorhombic (r/, rj, 0) and inverted hexagonal (r/, r/, rj) phases. Below, 
the phase digram stemming from Eqs. El El 

than that of the actual lamellar smectic phase, this decorrelation transition 
must be accompanied by vanishing of a physical quantity responsible for this 
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structural change. In other words the order parameter of the lamellar corre- 
lation transition can be defined. Group-subgroup relationship between these 
two lamellar phases permits one to show that the physical quantity vanishing 
through decorrelation transition is the macroscopic transverse polarization 
vector wave with the wave vector parallel to the smectic layer and the po- 
larization directed perpendicular to the single lamella (ID transverse vector 
wave). Condensation of this wave leads directly from the smectic- like decor- 
related structure to a correlated crystal lamellar phase of the rji = {rj, 0, 0) 
type, with periodic distribution of transverse polarization in columns in each 
layer. The order parameter of the transition is given by the first harmonic 
of the wave i.e. u^{k)ezexp{iky) + u^{—k)ezexp{—iky), where u^{k) is the 
amplitude of this harmonic and k is the wave vector perpendicular to the 
direction of the columns within a lamella, with the magnitude k = where 
a is the distance betweem neighboring columns within a single lamella. 

These polarization waves exist localy for each lamella, but the distribution 
of the origins of these waves in each lamella is random for the whole sample 
and consequently the average macroscopic value of the wave is zero in the 
decorrelated phase (Fig. 1). Similar scenario but in a different context was 
first proposed by P. Bak ^T]. The order-disorder character of the transition 
and the wave nature of the order parameter were clarified. It was also shown 
that the wave vector of the macroscopic order parameter for the transitions 
of this type is equal to the wave vector of the local waves (wave vector of 
an atomic chain in ^T] or of the polar vector wave in an individual lamella 
in our case). The transition from the decorrelated smectic phase into the 
correlated rji = {rj, 0, 0) type phase is described trivially via the Lifshitz free 
energy (T^j J-^ = ai\u''{k)\'^ + a2\u^{k)\'^ + . . . and is thus of second order. 

In this correlated-decorrelated transition scenario the crystalline lamellar 
phase is a common subphase of both the hexagonal crystalline phase as well as 
the smectic lamellar phase. The hexagonal crystal — > crystal lamellar phase 
transition is driven by condensation of the 2D periodic transverse polarization 
field, whereas the smectic lamellar — > crystal lamellar phase transition is 
driven by the condensation of the ID polar vector wave and is characterised 
by the correlation between the transverse polarization vector distributions in 
each individual lamella. 

Let us finally discuss the influence of chirality on the state of NCP 
columns. Each NCP particle being chiral their aggregation results in a chiral 
column with the rotational symmetry of a cylinder without neither mirror 
planes nor inversion, and the periodicity approximately equal to the indi- 
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vidual NCP height. According to the Curie principle ordering of the polar 
vector Pj along a chiral column leads to a simultaneous ordering of the axial 
vector tj, since both Pj and tj span the same representation of the column 
symmetry group. In a chiral periodic aggregate the axial vector has an evi- 
dent physical interpretation. Non-zero average value of tj corresponds to the 
correlated tilt of individual NCPs with respect to the column axis (Fig. 1). 
In the 2D hexagonal phase the NCP columns carry zero average polar vector, 
then in a column the tilt of the particles is also zero : the axes of individual 
particles (usually called S-axes) are parallel to the column axis. On the con- 
trary, in the lamellar phase the ordering of the polar dyadic axes induces a 
tilt of the NCPs S-axes. This simple consequence of the NCP chirality can 
explain peculiarities of the SAXS spectrum observed in jH]. 

One of the authors (R.P.) would like to thank F. Livolant and her collabo- 
rators for their hospitality while being a Chercheur Associe fellow of C.N.R.S 
at Orsay in 2002. 
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